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Recently, Mitsuhashi et al have observed superconductivity with transition temperature up to 
18 K in potassium doped picene (C22H14), a polycyclic aromatic hydrocarbon compound [Nature 
464 (2010) 76]. Theoretical analysis indicate the importance of electron-phonon coupling in the su- 
perconducting mechanisms of these systems, with different emphasis on inter- and intra-molecular 
vibrations, depending on the approximations used. Here we present a combined experimental and 
ab-initio study of the Raman and infrared spectrum of undoped solid picene, which allows us to 
unanbiguously assign the vibrational modes. This combined study enables the identification of the 
modes which couple strongly to electrons and hence can play an important role in the supercon- 
ducting properties of the doped samples. 
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Recently, superconductivity has been observed in potassium doped picene, with transition temperatures Tc of 7-18 
K^. This was the first report on " high- Tc" superconductivity (SC) in an aromatic compound. Picene (C22H14) is an 
ahernant polycychc aromatic hydrocarbon (PAH), i.e. a planar aromatic molecule, formed by juxtaposed benzene 
rings. Picene comprises five rings, arranged in a zig-zag fashion. After this initial report, SC was also reported in 
phenanthrene (three rings)^, dibenzpentacene (five + two rings)'^ and coronene (six rings)^ upon doping with alkali 
or alkalki-earth metals. SC was then also observed in picene and phenanthrene doped with rare-earths'*^^, indicating 
that the PAHs most likely form a completely new class of superconductors, with highly tunable electronic properties. 

Due to their 2-dimensional nature, PAHs are naturally related to other layered superconductors, such as MgB2^, 
or the newly-proposed graphane'*. However, in contrast to all these systems, which are conventional electron-phonon 
(ep) superconductors, the nature and origin of SC in PAH are still debated. Theoretical studies have focused either 
on electron-electron (ee)^'**^ or ep interaction**"*"*, which are both sizable in these 7r-bonded systems. However, 
estimates of strength and spectral distribution of the ep interaction differ depending on the approximation used for 
doping. The main question is to estabilish how much the ep coupling in an isolated molecule picture is representative 
of coupling in the actual doped compound**"*''. A definite answer on this point could only come from a detailed 
experimental comparison of the phonon spectra and lifetimes of the doped and undoped samples. This has not been 
possible so far, because the vibrational spectra are extremely rich, and the samples are poorly characterized. For the 
doped superconducting samples in ref.*, not only the detailed crystal structure, but also the precise doping level are 
unkwown. For pure picene. X-ray diffraction profiles have been refined several years ago *^, but, owing to the lack 
of high-purity samples, there were no serious attempts to a complete characterization of the vibrational spectrum, 
except for a recent study of the optical phonons**', which appeared in the literature while we were finalizing the 
present manuscript. 

In this paper, we present a combined experimental and theoretical study*^'*^ of the vibrational spectrum of crys- 
talline picene, based on accurate Raman and infrared measurements at room temperature on high quality samples. 
Density functional perturbation theory (DFPT*'') calculations were also performed, which allowed us to obtain a de- 
tailed assignment of vibrational modes. This study is propaedeutical to further studies on the high-pressure behaviour 
of doped and undoped samples, which could lead to significant insights into the importance of intermolecular coupling 
of vibrational modes and electrons in hydrocarbons, and on the reliability of DFT in organic compounds^*^ . Besides 
superconductivity, our study is of interest also for potential applications in organic electronics^'^*. 

Solid picene was prepared by a new optimized synthesis route which permits to obtain bulk quantities (hundreds 
of milligrams in one shot) of pure polycrystalline picene powder^^. Refinements of the X-ray diffraction spectrum of 
our sample yielded results identical to those previously reported*^. We performed Raman measurements by means of 
a confocal-microscope Raman spectrometer, equipped with a He-Ne laser (A=632.8 nm), a 1800 lines-mm^* grating 
monochromator, and a charge-coupled-device (CCD) detector. Raman measurements were carried out on the powder 
samples in the backscattering geometry, and a notch filter was used to reject the elastic contribution. Using a few 
micron-sized laser spot on the sample, accurate Raman spectra were collected with a spectral resolution better than 
3 cm^*. A typical Raman spectrum is reported in Fig. 1(a) over the frequency range 200-1700 cm^*. Note that 
intensity and width of the highest and sharpest peaks may be affected by the instrumental spectral resolution. For 
infrared (IR) measurements, we employed a thin pellet (of thickness around 50 micron) prepared out of pure picene 
powder. IR trasmittance measurements were carried out at the SISSI infrared beamline of the ELETTRA synchrotron 
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FIG. 1: (Color online) Comparison of the experimental Raman spectrum of crystalline picene (a) with the DFT spectrum 
of crystal (b) and relaxed molecule (c). For the molecule, the spectral components of different symmetry (ai, fei, 02, 62) are 
shown separately (see text for details). The computed DFT stick spectra are convoluted with a Lorenzian profile with 5 cm^^ 
linewidth to ease comparison with experiment. 



(Trieste, Italy) Using a Bruker IFS 66v interferometer equipped with both mylar beam splitters coupled to a liquid 
He bolometer for the far-IR region and a KBr beam splitter coupled to a mercury-cadmium-telluride (MCT) detector 
for the mid-IR region, accurate IR spectra with a spectral resolution close to 2 cm^^ were collected. A typical 
IR spectrum is reported in Fig. 2(a) over the frequency range 200-1700 cm~^. We notice that the most intense 
absorption lines result from vanishingly small transmitted intensities. These lines can therefore be affected by rather 
large uncertainties. Nevertheless, a large number of low intensity absorption peaks are determined with a good 
accuracy. 

The Raman and IR spectra of picene in its crystalline and molecular forms have been calculated using DFT^^'^'*'^''. 
The electron-ion interaction was described using LDA norm-conserving pseudopotentials^® ; electronic wavefunctions 
were expanded on a plane-wave basis, with a 80 Ryd cutoff. A 2"^ uniform grid for k-space integration in the solid 
allowed an accuracy of ± 5 cm^^ on the calculated frequencies. For the crystal, we employed the experimental 
unit cell, and relaxed atomic positions. The resulting Raman and IR spectra are reported in Figs. 1(b) and 2(b), 
respectively. For the molecule, we employed periodical repetitions of a 24 x 10 x 5 A supercell, containing a single 
picene molecule, with the z axis orthogonal to the ab plane, and the long molecular axis along the x axis. The cut-off 
energy was rescaled accordingly. We performed calculations both for an ideal geometry, and for a relaxed geometry 
in which we optimized the internal coordinates subject to the constraint of symmetry. The resulting Raman and IR 
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FIG. 2: (Color online) Same as Fig. 1, for the IR spectra. 



spectra are reported in Figs. 1(c) and 2(c), respectively. Note that the results we report in these figures are for the 
relaxed molecular geometry, but the difference between the ideal and relaxed molecule spectra is never larger than ± 
10 cm~^. 

In comparing the measured and the calculated spectra, it is convienent to recall that the calculated partial phonon 
density of states (not shown) indicates that modes up to 1000 cm^^ involve both out-of-plane and in-plane vibrations 
of carbon and hydrogen; above 1000 cm~^, the vibrations are only in-plane, with high-energy vibrations above 3000 
cm^^ involving only H atoms. The unit cell of crystalline picene is monoclinic (space group P2i) and contains 
two C22H14 molecules, arranged in a herringbone fashion and inclining slightly to the ab plane. This yields 216 
vibrational modes, many of which are both IR and Raman active. The comparison of the upper and middle panels 
of the two figures gives an idea of the accuracy of DFT calculations: DFT captures many of the main features of 
the experimental spectra, including the positions, intensities and splitting of the main peaks, with a few exceptions. 
The lower panel of the two figures shows the equivalent spectra, calculated for the picene molecule; the symbols show 
the total intensities, while colored lines refer to the four possible symmetries for the vibrational modes. The picene 
molecule has C2V symmetry; of the 102 optical modes, those with ai (35), 62 (16), 61 (34) symmetry are both Raman 
and IR active, while the 02 (17) are only Raman-active^^. A full list of peak frequency, symmetry, optical activity and 
the normalized intensity of the vibrational modes is given in the supplementary information file. It is worth noting 
that, for an isolated molecule, the Raman reponse is dominated by the (totally symmetric) ai modes (Fig. 1(c)), the 
IR reponse by the b2 modes (Fig. 2(c)). 

From a comparison between the calculated crystal and molecule spectra (see Figs. 1 and 2), it can be observed that 
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FIG. 3: (Color online) Comparison between the measured IR and Raman spectra (middle panel) and DFT Raman (upper 
panel) and IR (lower panel) spectra in the range 200 - 600 cm~^. The computed DFT stick spectra are convoluted with a 
Lorenzian profile with 10 cm^^ linewidth to ease comparison with experiment. Similar comparisons at higher frequencies are 
given in the supplementary information file. In the IR DFT spectrum of the molecule (lower panel), the 62 peak at 523 cm~^ 
marked by * splits into A and B components (see text) in the measured spectrum (middle panel). 



the molecular modes dominate the optical response. Picene is thus a typical molecular crystal, in which the optical 
response is determined by the selection rules due to the point group of the molecule, on which the site group of the 
crystal acts as a perturbation^^. 

Given the overall agreement between the spectra calculated for the molecule and the crystal, we have used the 
molecular symmetry assignments to label the vibrational modes observed in the solid. Fig. 3 shows a detailed 
comparison between the experimental (Raman and IR) and the calculated (molecule) spectra in the energy range 
between 200 and 600 cm^^. Similar comparisons at higher frequencies are given in the supplementary information 
file. In the same file we also report results obtained in the frequency range of the C-H stretching (3000-3250 cm^^), 
which is not of direct interest in the present work. Exploiting this kind of plots and of comparisons, we have been 
able to make a reliable assignments of most of the observed vibrational peaks. 

In Table 1, we provide the list of the most intense peaks which can be assigned without ambiguity. Our results 
extend a recent study^*^, which contains a partial assignment of the vibrational modes of picene, based on optical 
measurements and calculations'^^. Note that in correspondence to some of the DFT lines computed for the molecule, 
a doublet is observed; this effect can be attributed to the Davydov splitting of two modes corresponding to the same 
molecular vibration^®. For example, in Fig. 3 (see lower panel) the 62 peak at 523 cm^^ marked by * splits into two 
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TABLE I: Assignment of the most intense vibrational modes observed. For each mode, the frequency, symmetry and optical 
activity calculated in DFT are reported together with the corresponding IR and/or Raman frequency. All frequencies are in 
cm~^. The ai modes expected to be most involved in the ep coupling^'^ are marked by *. A complete list of frequencies and 
intensities of all the modes calculated in DFT is provided in the supplementary information file. 
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FIG. 4: (Color online) Distribution of the electron-phonon coupling at the F point, calculated in the rigid-band approximation 
in Ref.^^ and partial a\ phonon density of states. For details of the calculation see'^'\ 



components in the measured spectrum (see middle panel). In the table, two values corresponding to the same DFT 
mode indicate the observed splitting. 

The precise assignment of the phonon frequencies is particularly important, in connection with the problem of 
ep coupling, which in PAHs is usually calculated taking into account only vibronic contributions (see for example 
Refs.^'^'^''). For picenc doped with mono-, bi- and tri-anions, as well as with mono-cations, the coupling is non-zero 
only for modes with ai symmetry. According to Ref.^'', in the case of bi-anion doping, the modes most involved in 
the ep coupling are at 253, 724, 1380, 1516 cm"^. The corresponding values we observe for pure picene are 262, 
730, 1377 and 1516 cm~^ (marked by * in Table 1). A recent calculation by Casula et al. for K3 picene^^ has 
shown that this picture is strongly modified by the presence of the potassium, which mediates a strong intermolecular 
interaction. This result is strongly dependent on the geometrical distribution and concentration of the dopants, and 
is not representative of pure picene. 

In Fig. 4 we compare the spectral distribution of the ep coupling A, calculated in Ref.^^ using the rigid band 
approximation, with the partial ai phonon density of states to estimate the effect of intermolecular interactions of 
pure crystalline picene'''^. The figure shows that the ep coupling response is indeed dominated by the ai modes, but 
there are additional contributions clearly coming from modes with different symmetries. Consistently with what we 
observed in the IR and Raman spectra, the figure shows that intermolecular interactions in undoped picene are quite 
small; this picture could be modified by doping and/or pressure, which could both induce a stronger intermolecular 
coupling, due to hybridization and/or increased overlap between the molecules. 

In conclusion, we have presented high-quality optical data and first-principles calculations for the vibrational 
spectrum of solid picene. The comparison between experimental and theoretical results enabled us to completely 
characterize the spectrum, and assign most of the observed vibrational modes. Comparing the theoretical calculations 
for Raman, IR and electron-phonon spectral function, we have identified picene as solid, with a small intermolecular 
interaction. The extensive characterization obtained in this study is fundamental for further studies on the ef- 
fect of pressure and/or doping on bonding and intermolecular interactions in doped polycyclic aromatic hydrocarbons. 
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DFT phonon modes in the frequency range 200-1700 cm^. For each mode, calculated DFT frequency, 
symmetry, optical activity, and intensity of the IR and/or Raman frequency are reported. Intensities are 

normalized to a maximum value of 100.00. 
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FIG. 5: Comparison between the measured IR and Raman spectra (middle panel), DFT Raman (upper panel), and IR (lower 
panel) spectra in the range 600 - 1100 cm~^. The computed DFT stick spectra are convoluted with a Lorenzian profile with 
10 cm^^ linewidth to ease comparison with experiment. 
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FIG. 6: Comparison between the measured IR and Raman spectra (middle panel), DFT Raman (upper panel), and IR (lower 
panel) spectra in the range 1100 - 1700 cm^^. The computed DFT stick spectra are convoluted with a Lorenzian profile with 
10 cm~^ linewidth to ease comparison with experiment. 
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FIG. 7: Comparison between the measured IR and Raman spectra (middle panel), DFT Raman (upper panel), and IR (lower 
panel) spectra in the frequency range of the C-H stretching (3000-3250 cm~^) which is not of direct interest in the present 
paper. Notice that the measured Raman and IR spectra have been normalized to the same maximum. Like in other frequency 
ranges, the computed DFT stick spectra are convoluted with a Lorenzian profile with 10 cm^^ linewidth to ease comparison 
with experiment. 



